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Abstract. We probe the validity of Crooks' fluctuation relation on the fluctuating 
lattice-Boltzmann model (FLBM) , a highly simplified lattice model for a thermal ideal 
gas. We drive the system between two thermodynamic equilibrium states and compute 
the distribution of the work performed. By comparing the distributions of the work 
performed during the forward driving and time reversed driving, we show that the 
' system satisfies Crooks' relation. The results of the numerical experiment suggest that 

the temperature and the free energy of the system are well defined. 
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(Ti ■ 1. Crooks Fluctuation Relation 

^ I The second principle of thermodynamics states that on average one cannot extract 
O I work from a thermodynamical system coupled with only one heat bath during a cyclic 
^ ' process. To illustrate this, consider a thermodynamical system that we can drive with a 
parameter A, in contact with a heat bath at constant temperature. By switching A from 

to 1, we can drive the system from one equilibrium state to another, and by doing this, 
^ I we will perform some work Wf. Consider the time reversed or backward experiment: 

By switching A from 1 down to 0, we drive the system back to its initial state 1§| and we 
extract some work — W]^. The second law of thermodynamics states that on average we 
extract less work during the backward experiment than we perform during the forward 
experiment: 

(Wi) > -{WQ. (1) 

1 Present address: Carl von Ossietzky Universitat Oldenburg, Institut fiir Physik, 26111 Oldenburg, 
Germany 

§ By state we mean macroscopic equilibrium state specified by the value of A, the temperature of the 
heat bath and other external constraints. In fact, we assume that, given a set of external constraints, 
the system will always relax towards thermodynamic equilibrium. 
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The equahty in the above inequahty holds if and only if the forward and the backward 
process are performed quasi- statically, i.e., infinitely slowly such that the system always 
remains in equilibrium with the heat bath. In that case, the work performed during 
the forward process or extracted during the backward process is called the reversible 
work and is equal to the difference in free energy between the final and the initial states 

Wt = AF = Fi- Fq. (2) 

The quantity = Wf — AF (resp. = AF + W\y) is then the dissipated work during 
the forward (resp. backward) process. 

Crooks' fluctuation relation [1] gives some quantitative information on the 
probability to dissipate a certain amount of work during a non quasi-static process. 
Consider a particular forward protocol, specified by the time dependence Af(t) of the 
driving parameter A. Now also consider the time reversed or backward protocol obtained 
from the former simply by inverting the time dependence of A: )^i){t) = Af(ts — t), 
where ts is the switching time. Crooks' fluctuation relation states that the ratio of 
the probability Pf(W^^) of dissipating a certain amount of work VF^ using the forward 
protocol and the probability P-\^{—Wf^) of dissipating the opposite amount —W^ using 
the time reversed protocol is given by: 

^fC^d) -expfi%V (3) 



where ks is Boltzmann's constant and T is the temperature of the heat bath. Note that 
the dissipated work may be negative. We expect it to be positive only on average. Given 
the definition of W^, the ratio on the left hand side of (|3]) is the same as the ratio of the 
probability to perform a certain amount of work W using the direct protocol over the 
probability to perform the opposite amount of work using the time reversed protocol. 
Crooks' relation then reads: 

A sufficient condition for the validity of this relation is that the system satisfies detailed 
balance for fixed A [T]: 

P{A ^ B) _ e-^^B 

P{A ^B)~ e-P'^A ' ^ ' 

where P{A — )■ B) is the probability that the system in microstate A with energy E\ 
makes a transition to microstate B with energy and /3 = l/fc^T is the mverse 
temperature of the heat bath. This condition is very restrictive and can be relaxed to 
the less restrictive condition of balance, see [?] for more details. 

Since the discovery of the Jarzynski equality in 1997 [2], fluctuation theorems have 
been intensely studied experimentally [H H] and theoretically [5] . The purpose of this 
work is to probe the validity of Crooks' fluctuation relation on a very simple, but realistic 
lattice gas model. 

The model chosen for the numerical experiment is the fluctuating lattice-Boltzmann 
model (FLBM). It constitutes one of the simplest lattice models for an ideal thermal 
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gas. Classical lattice-Boltzmann models are a powerful tool for the simulation of 
hydro dynamic systems [6]. However, since they usually simulate systems on macroscopic 
scales, they do not include thermal fluctuations. The FLBM was developed in order to 
simulate the solvent in soft matter applications, where thermal fluctuations are not 
negligible [3 IE]. In this work, we simulate such an ideal gas subjected to a potential 
field per unit mass 4>\{r) depending on a certain control parameter A, and we measure 
the work performed on the system by switching A from to 1 and from 1 down to in 
a finite time ts- 

Probing Crooks' relation on the FLBM is a way to check the thermodynamic 
consistency of this model. In fact, the FLBM was developed in order to correctly 
simulate equilibrium fluctuations, but has no clearly defined total energy or free energy. 
By probing Crooks' relation, we mean checking that the ratio Pf{W) / P]^{—W) is an 
exponential function of W and that this function is independent of the specific protocol, 
that is the time dependence of the driving parameter A. We can then check that the 
exponent /3 appearing in this exponential function is the inverse temperature l/ksT 
of the the heat bath. The main result of the present paper is that this model indeed 
satisfies Crooks' relation (jl]) in a thermodynamically consistent manner. After having 
briefly introduced the model, we will present the numerical experiment. Finally, we will 
focus on the distribution of the work performed during the switching and show that 
for this particular case, the system satisfies Crooks' fluctuation relation. Moreover, the 
distribution of work turns out to be Gaussian, which permits analytical simplifications 
of (H. 

2. The fluctuating lattice-BoItzmann model 

Our purpose here is not to give a detailed description of the fluctuating-lattice 
Boltzmann model (FLBM). We will simply highlight the main features of the model 
we are interested in and we refer the interested reader to the literature for more details. 
A good introduction to lattice-Boltzmann models can be found in [6] and the particular 
model used here is described in detail in [TJ IE]- The FLBM is a highly simplified gas 
model. Its dynamics takes place in discrete time steps and on a square lattice. At each 
time step, the particles move from one node of the lattice to another according to a 
discrete set of velocities and then they collide with the particles sitting on the same 
node. The time step and the lattice spacing are taken to be unity defining the time and 
length units. The mass unit is arbitrary, but is the same throughout this document. 
The mass of one particle per unit volume is denoted by yU such that if p is the mass 
contained in one unit volume, then is the number of particles contained in this unit 
volume. Throughout this document, all the physical magnitudes are expressed in this 
unit system. 

The model used for this work is the D2Q9 model. It consists of 9 velocities on 
a two dimensional square lattice, (see figJI]). Let us denote by {cj,z = 0, ...,8} the 
set of velocities. The dynamical variables of this model are the set of populations 
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Figure 1. The two dimensional nine velocities set of the D2Q9 model: Cq is the zero 
velocity, Ci, . . . , C4 connect to the nearest neighbors and C5, . . . , Cg connect to the next 
nearest neighbors. 



n = {ui, i = 0, . . . , 8} of the different velocities. The occupation number nj(r, t) is the 
mass density on node r = {x, y) moving along velocity Cj at time t, where {x, y) are 
the Cartesian components of r. The mass density p and the mass current density j 
are given by the zeroth and first moments of the velocity set with respect to the set of 
populations: 

8 

Pir,t) = Y.Mr,t), (6) 

2=0 

and 

8 

j{r,t) = Y.n,{r,t)ci. (7) 

i=0 

The local velocity field v is then given by v{r,t) = j{r,t)/ p{r,t). 

The dynamics of the system takes place in two steps: propagation and collision. 
During the propagation step, the particles simply move from one node to the other 
according to their velocity. During the collision step, all the particles sitting on 
the same node interact with one another exchanging momentum according to some 
basic conservation rules. The algorithm can be summarized by the lattice-Boltzmann 
equation: 

mir + Ci,t + 1) = n*{r,t) = ni{r,t) + Ai{n) + Fi. (8) 

Here, n*{r,t) is the occupation of velocity Cj just after the collision step and {Aj} is 
the collision operator. The last term Fi simulates the action of a force per unit volume 
F applied on the system. The collision operator and the term due to the force will be 
presented just after a discussion of the equilibrium fluctuations. 

2.1. Equilibrium fluctuations 

As mentioned before, in classical lattice-Boltzmann simulations, the collision operator 
is deterministic and simply lets the velocity distribution {rij} relax towards its local 
equilibrium value {rif^} [8j. The local equilibrium populations {rif*^} are functions of 
the locally conserved quantities p and v that we will specify later. In addition to this 
relaxation, the FLBM allows for small fluctuations around local equilibrium. These 
small fluctuations are assumed to be Gaussian in first approximation. Taking mass and 
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momentum conservation into account, the equilibrium probability distribution of the 
velocity occupations is [7]: 

' i:fe-iME".-.iME"..-.-).(o) 



P{{ni}\p,v) oc exp 



v 




where 5{x) = 1 for a; = and otherwise. The equilibrium fluctuations are controlled 
by /i: 

{nj) - {n,f = fia^^p, (10) 

The idea behind this relation is that the fluctuations in the number nj//i of particles 
having velocity Cj at a given node is proportional to the total number of particles p/fi 
on that node. The weights a'^' serve to restore isotropy in the large scale limit. They 
must sum up to unity, = 1; and must be compatible with the symmetries of the 

lattice [TJ |8] and in particular, they do not depend on the direction of the velocity Cj 
but only on its magnitude. For the D2Q9 lattice used here, their values are: a° = 4/9, 
= 1/9, = 1/36 (see p.69). 
The particle mass per unit volume p is our fluctuation parameter. It is linked to 
the resolution of the simulation. The limit p corresponds to the thermodynamic 
limit: decreasing p means to describe the system on a more coarse grained scale so that 
the number of particles we describe goes to infinity and the amplitude of the thermal 
fluctuations of the macroscopic observables decreases to zero. On the other hand, if p 
is of order one, then the fluctuations of the mass density are of the same order as the 
mean. 

The temperature T of the system is proportionnal to the fluctuation parameter p: 
ksT = pel (11) 



Here, cs = \Jdp/dp is the isothermal speed of sound. Equation (|TT1) is a consequence of 
the equation of state of ideal gas[Jj] 

P = -ksT, (12) 
p 

where p is the hydrodynamic pressure of the gas. The speed of sound cs is an intrinsic 
quantity of the set of velocities : it is the maximum speed at which a signal can propagate 
through the system. It is of the order of |cj|. For the D2Q9 model, its value in lattice 
units is Cs = l/V^ ([6] p.69). 

II Note that (|12|) and are inhomogeneous. Indeed, /i is a mass density: ^ — rrip/b'^, where nip 
is the mass of one particles and b is the lattice spacing. But, by considering b ~ 1, we write /i ~ rup 
for simplicity of the notation. Moreover, we refer to /i as the fluctuation parameter and not as the 
particle mass or as the temperature to keep in mind that it is a coarse graining parameter: Changing 
fi definitely changes the temperature of the system according to ([TT|). but it also changes the number 
of particles through rrip. For a discussion about fj,, see [9]. 
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2.2. The collision operator 

The collision operator {Aj} computes the changes in the populations due to the 
collisions: 

A,(n) = <(r,t)-ni(r,t), (13) 

where ni{r,t) is the precollisionnal value of the occupation of velocity Cj at node r and 
time t and n*[r,t) is its postcoUisionnal value. The role of the collision operator is to 
thermalize the system. It operates as follows: The velocity distribution {rij} at each 
node is linearly relaxed towards its local equilibrium value {nf*^} and this relaxation is 
balanced by a suitable thermal noise. The equilibrium populations {^f*^} depend only 
on the local values of p and j, and are such that mass and momentum are conserved: 



Enl'' = P (14) 

i=0 



Y^n^%i=j, (15) 

i=0 

and such that the equilibrium stress is correctly given by the Euler stress 

8 

n^%iaCii3 = pSal3 + pVaVp, (16) 

i=0 

where the hydrodynamic pressure is given by the ideal gas equation of state (fT2l) 
p = -ksT = pcg, and a and /3 are the Cartesian coordinates. Those conditions uniquely 
determine the local equilibrium populations. They read [7]: 

The postcoUisional state n* is obtained from the precoUisional state Ui by the 
collision rule 

n* - = (rii - n^^^ + Ti, (18) 

where 7 is a relaxation parameter satisfying I7I < 1 linked to the relaxation time of the 
velocity distribution, and rj is a random number corresponding to the fluctuating part 
of the collision process. 

The conservation of mass and momentum imposes 

E^. = 0, (19) 

i=0 

and 

8 

Y.r,c, = (20) 

i=0 

% This collision rule is the same as in [71 15] in the case of one single relaxation time. 
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such that, in the end, there are only 6 independent random numbers corresponding 
to the 9 degrees of freedom minus 3 constraints. The rjS are drawn from a Gaussian 
distribution with zero mean, and their variance is 

(r^) = /.a^^py,, (21) 

with if = y/1 — 7^ so that the variance of is correctly given by (ITO!) . It can be shown 
[7] that the collision rule ( ITSl) satisfies detailed balance with ([9]) as invariant distribution. 



2. 3. Applying a body force to the system 

In order to drive the system, we apply a force f{r,t) per unit mass on the system. The 
force density applied on one lattice site is then F = pf. The effect of the force is to 
increase the momentum density j by 6tpF, at each time step, where 5t = 1 is the time 
step. This is done by adding a deterministic term Fi = pfi to the lattice-Boltzmann 
equation ([8]). The fi must satisfy mass and momentum conservation, namely 



E = 0, (22) 



i=0 



and 



Y.f^<'^ = f- (23) 

As mentioned in [8], the definition of the local velocity v is no longer unique due 
to the discretization of the time. Its value at each time step is included between its 
precollisional value J2i ^i^i / J2i and its post-collisional value + /• As 

in [8j, we define the local velocity as 

« = -^n,c, + l/, (24) 
P i=0 ^ 

which is the mean value between the precoUisionnal and the postcoUisonnal values. 

In [8], the expression of Fi consistent with hydrodynamics is derived. Following 
from that, the expression for the fi is 

1 + 7 



f^ 



+S 




■ /• (25) 



3. Numerical experiments and results 

We drive our system form one equilibrium state denoted by to another denoted by 
1 according to a particular protocol X(t), A(0) = and A(ts) = 1, where is the 
switching time. The backward process is obtained by starting in the equilibrium state 1 
and tuning A down to 0. For each run, we record the work W performed. By repeating 
the forward and the backward experiments many times, we get the distribution Pf{W) 
of the work performed during the forward process and the distribution P]j{—W) of the 
work extracted during the backward process. 
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Figure 2. Potential applied on the system as a function of x expressed in units of 
lattice spacing. The forward process is started with A = 0. At each time step, A is 
incremented by a small amount until A = 1. 



3.1. Direct and time-reversed protocol 

The macroscopic state of the system is controlled by applying a force field per unit 
mass fx(t){f) = — V0A(t)(T') derived from a parameter dependent potential (f>x{r) on the 
system. The potential energy of the system is 

r,(t)=^p(r,t)0,(,)(r). (26) 



The work performed on the system during one run is 



W 



*s 6W dX 

1 

5\ dt 



-dt = ± 



6W 



dA, 



(27) 



where the plus sign holds for the forward process and the minus sign for the backward 
process. The work performed at each time step SW is the variation of potential energy 
of the system when A is varied by a small amount 6X, while no heat is exchanged with 
the heat bath: 



SW = Vx+5x -Vx = i<f>x+5x{r) - Mr)) 

r 

The potential used in this work is (see figj2]): 



XA 



cos ( 2ttj ] + 1 



(28) 



(29) 



where / is the length of the system and A is the amplitude of the potential. Thus, the 
force applied to one unit of mass is 



fx 



27r 



V0A = XA— sin 27r— e^. 



X 



I 



For the direct protocol X{t) = t/ts and for the time- reversed protocol X(t) = 1 
The work 6W performed during one time step is then 



6W{t) = ASXj2p{r,t) 



X 



cos I 2nj I + 1 



(30) 
t/ts. 

(31) 



where 6X 



dX 
dt 



±l/ts is positive for the forward process and negative for the 



corresponding backward process. 




Figure 3. Distribution of the work for three values of the switching rate and for 
fi = 10. The thin hues correspond to the work performed during the forward process 
and the thick hues to the work extracted during the backward process. The work is 
given in lattice units. 



The simulations were performed on a / x L = 100 x 10 lattice with an average 
density of 1000 and periodic boundary conditions. The relaxation parameter 7 was set 
to 0.9 and the amplitude A of the potential to 0.01. The parameters of the model are 
/i, controlling the amplitude of the fluctuations, and the switching rate, 6X, controlling 
how far we are from a quasi-static process. 

3.2. Fluctuations of the work performed 

An example of the distribution of the work performed during the forward process and the 
work extracted during the backward process can be seen on figure [31 These simulations 
were performed with /x = 10 for various values of SX. According to Crooks' relation (jl]), 
the reversible work AF is the value of W for which Pf{W) = p\j(—W). As can be seen 
in figure [3l this value is independent of the switching rate 6\, as expected. Moreover, we 
can already see that, as the switching rate decreases, the mean work converges towards 
the reversible work AF and the fluctuations of the work go to 0. In the limit 6X — ?■ 0, 
we expect the work distributions to be Dirac distributions peaked at AF, as predicted 
by equilibrium thermodynamics. 

3.2.1. Influence of the fluctuation parameter The fluctuation parameter /i controls the 
amplitude of equilibrium fluctuations of p according to (fTOl) . During the driving, the 
situation is not very different. To illustrate this, we plotted examples of the evolution of 
the work 6W/SX performed at every time step during the forward process as a function 
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Figure 4. Work performed at each time step as a function of A for three different 
reahzations of the forward process with different values of the fluctuation parameter 
/i. Simulations performed with 6\ — 10^*. 

of A for three different values of /i on figure HI According to equation ( j27l) . the total 
work W performed during the switching is the integral of this function. 

As we can see on this figure, the amplitudes of the fluctuations of 5W are clearly 
controlled by ^. For /i = 0, the evolution of the system is completely deterministic and 
5W does not fluctuate, and the bigger n is, the larger are the fluctuations. On figure 
O we plotted the distribution P{{W) of the total work performed during the forward 
process for various values of /i. On this figure, we can see that /i controls the fluctuations 
of the total work performed during the forward process without influencing its mean 
value. In fact, the variance o"^ of the distribution of the work is proportional to /i (see 
top left of fig. [5]). Moreover, as we will see later, the variance of the work performed 
during the backward process is the same as for the corresponding forward process. 

3.2.2. Influence of the switching rate To study the influence of the switching rate on 
the mean work performed during the process, let us first consider the case without 
fluctuations (/i = 0). Figure [6] shows dWjdX as a function of A for = for various 
switching rates. In that case, the work performed during the forward process is always 
bigger than the work extracted during the backward process. The difference between 
both is linked to the work dissipated during the process and it goes to zero as the 
switching rate goes to zero, that is as we approach a quasi-static process. 

For a finite amplitude of the fluctuations (/i > 0), the total work performed during 
the experiment fluctuates less for a slow switching than for a fast switching (see figj3]). 
However, the amplitude of the fluctuations of ^ are independent of the switching 
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Figure 5. Distributions of the work performed during the forward process for 



(5 A = 10 ^ and for different values of ji. 
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Figure 6. Work performed during the forward process (thin lines), and extracted 
during the backward process (thick lines) at each time step as a function of the driving 
parameter A in the absence of fluctuations (/z = 0). 




Figure 7. Fluctuations of tlie work performed at eaclr time step for two realizations 
of the experiment. Simulations performed with /i — 10. 



rate SX. This can be seen on figure [TJ where we plotted the fluctuations in the work 
performed at every time step as a function of A. The mean work (SW) was 

obtained from a simulation with /i = 0. The difference in the fluctuations of the total 
work W comes from the fact that for a slow switching, the integral fl27|) is computed 
over a longer time than for a fast switching, such that the fluctuations in ^ are partly 
integrated out. For a quasi-static process, the total work W does not fluctuate because 
the integral fl271) is computed over an inflnite time such that the fluctuations in ^ are 
completely integrated out. 



3. 3. Verifying Crooks ' relation 

Since we perform a flnite amount of realizations of the experiment, we can probe 
Crooks' fluctuation theorem empirically only for values of W where Pf{W) and Pi){—W) 
signiflcantly overlap. On flgure [H we plotted the ratio pf(W)/p]^{—W) for fi = 10 for 
various values of the switching rate. The line represents what we would expect from 
Crooks' relation, namely exp [/3{W — AF)]. The temperature ksT = = /ic| is given 
by (fTTj) and is 0.3 for fi = 10 and the reversible work AF is the value of W for which 
Pf{W) = p\^{—W) and can be obtained from flgure [3l When both Pf(W) and P\^{—W) 
are sufficiently large, their ratio satisfies Crooks' relation @. 

As can be seen on figure [9], the distribution of the work is Gaussian. As pointed 
out in [3], in that case we can simplify Crooks' relation. If we write the distribution of 
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Figure 8. The ratio p^{W) / p^{—W) as a function of for ^ = 10. 



the work performed during the forward and the backward process as a Gaussian 

2n 



1 



: exp 



W -W 



f,b 



2-f,b 



(32) 



then, the ratio appearing in Crooks' relation becomes 

J!fW . ^exp f (-L 1 ) , (l^W, . ^W,) H' + ^ - 1|) . (33, 

and Crooks relation (jl]) can be valid only if the variances for the forward process and 
the backward process are equal: = = a"^, so that the term in W"^ vanishes in the 
exponent and that the factor in front of the exponential is 1. By identification, we can 
express the inverse temperature /3 and the reversible work AF in terms of the means 
and the variance of the work distributions: 



and 



(3 



AF 



^f+(-^b) 



(34) 



(35) 



These relations enable us to probe the validity of Crooks' relation in the Gaussian case 
even if the distributions for the forward process and the backward process don't overlap. 

The dissipated work (W^ = Wf — AF for the forward process and = W]^ + AF 
for the backward process) is then also Gaussian, and its mean is related to its variance 
by 

Wa = (36) 



Crooks ' fluctuation theorem for the FLBM 
1 I ^ ^ ^ — 



14 




le-06 I ^ ^ ^ ^ ^ ^ ^ 1 

-6 -4 -2 2 4 6 

Figure 9. Distribution of = W^/a — fia/2 for tlie forward process and the 
backward process. The data plotted come from simulations with = 10, and 
5\ = 10-4, 10-3, 5 • 10-3,^nd 10-2, and for 5\ = 10-^ and M = 3, 5, and 7. The line 
represents l/-\/27r exp(— i Wj). 



This relation can be seen as a generalization of the equilibrium fluctuation dissipation 
theorem. The mean dissipation is proportional to the fluctuation of the dissipation and 
the proportionality coefficient is the temperature, just as in the equilibrium fluctuation 
dissipation theorem. 

As a consequence of eq. ( l36ll . the quantity 

^ a 

is Gaussian with zero mean and unit variance, independently of the values of 6\ and /i. 
On figure [H we can see the distribution of for various values of /i and 5A, for the 
forward and for the backward process. As we can see on this figure, the distributions of 
coming from simulations performed with different values of the switching rate and of 
the fiuctuation parameter all collapse onto a Gaussian distribution with zero mean and 
unit variance From that, we can conclude that P{{W) and P\){—W) are Gaussian, that 
the system satisfies Crooks' relation (jlj) or equivalently ( l36l) . and that the temperature 
fc^T = and the reversible work AF are correctly given by Crooks' relation or 
equivalently by equations ( l3l|l and ( 135|) . 

4. Summary and outlook 



We have presented a basic numerical experiment which permits to probe the validity of 
Crooks' fiuctuation relation on the fiuctuating lattice-Boltzmann model. We have seen 
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that this model fulfills Crooks' relation and that in this particular situation, Crooks' 
relation is considerably simplified due to the Gaussian nature of the distribution of the 
work W performed during the switching. This simplification enables us to probe Crooks' 
relation, even in situations where Pf{W) and P\,{—W) do not overlap significantly. 

This experiment suggests that the FLBM is thermodynamically consistent. The 
reversible work AF depends only on the initial and final states, that is on A(0) and 
A(ts)- The mean dissipated work depends only on the switching rate 5\. The 
fiuctuations depend on the fiuctuation parameter yU and on the switching rate 5\. 
Even though this model has no a priori well defined energy or free energy, one could 
define the free energy (up to an additive constant) using fl35p . However, the link between 
this quantity and the microdynamics of the system is still to determine. 

The fiuctuating lattice-Boltzmann model is a very simple model for a thermal gas. 
It enables us to study Crooks' fiuctuation relation in a simple situation. Equation ( l27|l 
suggests that the interesting quantity is Its value during the process only depends 
on the instantaneous mass density distribution p(r', t). One could, for instance, study the 
fiuctuations in p{r,t) around its mean value during the switching and compare them to 
the equilibrium fluctuations. More generally, this model could help identify the sources 
of dissipation in a finite switching-time experiment or in a non-equilibrium steady state 
and the link between fiuctuation and dissipation in non-equilibrium situations. 
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